Competition between hidden order and antiferromagnetism in URu2Si2 under uniaxial 

stress studied by neutron scattering 
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We have performed elastic neutron scattering experiments under uniaxial stress a applied along 
the tetragonal [100], [110] and [001] directions for the heavy electron compound URu2Si2. We found 
that antiferromagnetic (AF) order with large moment is developed with o along the [100] and [110] 
directions. If the order is assumed to be homogeneous, the staggered ordered moment fi continu- 
ously increases from 0.02 iis/U (a = 0) to 0.22 /is/U (0.25 GPa). The rate of increase dfi /da is 
~ 1.0 ii_g/GPa, which is four times larger than that for the hydrostatic pressure (d/j, /dP ~ 0.25 
Hb /GPa). Above 0.25 GPa, \i shows a tendency to saturate, similar to the hydrostatic pressure 
behavior. For a\\ [001], fi shows only a slight increase to 0.028 \xb /U (a = 0.46 GPa) with a 
rate of ~ 0.02 tis/GPa, indicating that the development of the AF state highly depends on the 
direction of a. We have also found a clear hysteresis loop in the isothermal /x D (o") curve obtained for 
cr 1 1 [110] under the zero-stress-cooled condition at 1.4 K. This strongly suggests that the er-induced 
AF phase is metastable, and separated from the "hidden order" phase by a first-order phase tran- 
sition. We discuss these experimental results on the basis of crystalline strain effects and elastic 
energy calculations, and show that the c/a ratio plays a key role in the competition between these 
two phases. 

PACS numbers: 71.27.+a, 75.25.+Z, 75.30.Kz 75.30.Mb 



I. INTRODUCTION 



The nature of the phase transition at T a = 17.5 K in 
URu 2 Si2 (the ThCr 2 Si2-type, body-centered tetragonal 
structure is presently one of the most challenging 
issues in the heavy-electron physics. The elastic neu- 
tron scattering experiments^^ indicate that the simple 
type-I antiferromagnetic (AF) order develops below ~ 
T Q . However, the obtained staggered moment /i is ex- 
tremely small (~ 0.02 — 0.04 /zb/U), and incompatible 
with the large bulk anomalies such as the specific heat 
jump at T (AC/T ~ 300 mJ/K 2 mol). This incon- 
sistency has been puzzling many researchers for almost 
twenty years, i.e., whether the intrinsic order parameter 
is the tiny magnetic moment or some unidentified "hid- 
den" degree of freedom. The key to this issue has been 
recently obtained from the microscopic studies performed 
under hydrostatic pressure P. We found from the neu- 
tron scattering experiments that /z D is strongly enhanced 
by applying pressure from 0.017 /zb/U (P — 0) to 0.25 
/zb/U (P = 1.0 GPa)»££ In parallel, the 29 Si-NMR study 
revealed that the system is spatially separated into two 
differently ordered regions below T Q : one is AF with a 
large moment and the other is non-magneticAiS The AF 
volume fraction is found to increase with P, roughly in 
proportion to /z 2 (P) , while the magnitude of internal field 
is almost independent of P. This indicates that the ob- 



served enhancement of the AF Bragg-peak intensities is 
attributed to the increase of the AF volume fraction, and 
not of the local AF moment. Simple extrapolation yields 
the AF volume fraction at ambient pressure of about 1%, 
strongly suggesting that this is the true nature of the tiny 
magnetic moment. Consequently, the remaining 99% is 
considered to be occupied by the "hidden order" , which 
is responsible for the large bulk anomalies at T a . 

The major purpose of the present study is to inves- 
tigate how these two types of order correlate with each 
other. In order to find a relevant parameter, we here 
examine the effects of lattice distortion. So far various 
ideas for the hidden order parameters have been pro- 
posed, including valence transition^ uranium dimers^a 
unconventional spin density waves quadrupolar 
oj-denASiASiiLiSiASL and charge current order All of 
them involve a magnetic instability such that the dipolar 
order may be replaced with the majority hidden order. 
This switching is expected to be driven by lattice dis- 
tortion, since the proposed hidden order parameters are 
tightly coupled to the lattice system. It is thus interest- 
ing to investigate the competition between the two types 
of order by tuning the crystal distortion. 

A second purpose is to find the relationship between 
the two ordered states. The 29 Si-NMR results indicate 
that the AF order develops in parts of the crystal. How- 
ever, it is not clear whether it is inevitably induced 
through some coupling with the hidden order parame- 
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ter, or simply replaced with hidden order by a first order 
phase transition. In the latter case, hysteretic behav- 
ior can be expected in the pressure variations of the AF 
state. This point, however, was not checked in the pre- 
vious measurements,''' 8 where samples were always com- 
pressed at room temperature. 

For these purposes, we have performed elastic neutron 
scattering experiments on URu2Si2, by applying uniax- 
ial stress cr under both the stress-cooled and the zero- 
stress-cooled conditions. We have previously reported 
some experimental results obtained for weak a up to 0.46 
Qp a( 22,2£ j n j.jjg p resen ^ paper, we have extended the a 
range up to 0.61 GPa, and also investigated a Rh-doped 
system U(Ruo.99Rho.oi)2Si2- The collected results are 
discussed and interpreted in terms of a lattice distortion 
(or stress) model involving a distribution of the c/a ratio. 

II. EXPERIMENTAL PROCEDURE 

A single-crystalline sample URu2Si2 was grown by the 
Czochralski pulling method using a tri-arc furnace, and 
vacuum-annealed at 1000°C for a week. Three plates 
with three different bases of (001), (100) and (110) planes 
were cut from the crystal by means of spark erosion. The 
dimensions of the plates are approximately 25 mm 2 x 1 
mm. The uniaxial stress a was applied along the [001], 
[100] and [110] axes up to 0.61 GPa, by placing the sam- 
ples between Be-Cu piston cylinders mounted in a clamp- 
type pressure cell. This cell was used for measuring the 
temperature variations of the AF state down to 1.5 K 
under the stress-cooled condition, where the stress was 
changed at room temperature. 

The elastic neutron scattering experiments were per- 
formed by using the triple-axis spectrometer GPTAS 
(4G) located in the JRR-3M research reactor of Japan 
Atomic Energy Research Institute. The neutron momen- 
tum k = 2.660 A~ was chosen by using the (002) reflec- 
tion of pyrolytic graphite (PG) for both monochromating 
and analyzing the neutron beam. We used the combina- 
tion of 40'-80'-40'-80' collimators, together with two PG 
filters to eliminate the higher order reflections. The scans 
for the stress-cooled process were performed in the (hkO), 
(hOl) and (hhl) scattering planes for a || [001], [100] and 
[110], respectively. The AF Bragg reflections were ob- 
tained by the (100) scans for a\\ [001], the (100), (102) 
and (203) scans for cr 1 1 [100], and the (111) and (113) 
scans for a || [110]. 

For the measurements under the zero-stress-cooled 
condition, we used a constant-load uniaxial stress 
apparatus. 24 In this apparatus, the Be-Cu pistons in the 
pressure cell, which is attached to the bottom of the 4 He 
cryostat insert, is compressed by an oil-pressure device 
mounted on the top of the insert via a movable rod made 
of stainless-steel and tungsten carbide. The load is pre- 
cisely stabilized by controlling the oil pressure during the 
measurements. We first cooled the sample down to 1.4 
K without compression, and then applied the uniaxial 



stress along the [110] direction up to 0.4 GPa, keeping 
the sample at the same temperature. The scans for the 
zero-stress-cooled condition were performed in the (hhl) 
scattering plane. The AF Bragg reflections were ob- 
tained by the longitudinal scans at the (111) position. 
The experiments under the zero-stress-cooled condition 
(cr 1 1 [100]) were also performed on the Rh-doped alloy 
U(Ruo.ggRho.oi)2Si2, which was prepared in the same 
procedure as the pure compound. The (100) magnetic 
Bragg reflections were investigated by using longitudinal 
scans in the (hOl) scattering plane at 1.4 K. 



III. EXPERIMENTAL RESULTS 

A. Elastic neutron scattering under stress-cooled 
condition 

Figure^shows the a variations of the longitudinal and 
transverse scans at 1.5 K through the (100) magnetic 
peak for cr || [100] and || [001], and the longitudinal scans 
through the (111) peak for a || [110]. The instrumental 
background and the contamination of the higher-order 
nuclear reflections were carefully subtracted by using the 
data taken at 40 K. As stress is applied along the [100] 
direction, the (100) peak intensity markedly increases 
(Fig. da) and (b)). The (102) and (203) peaks also de- 
velop rapidly (not shown). The intensities of these three 
magnetic reflections divided by the polarization factor 
roughly follow the |Q| dependence of the magnetic form 
factor 25 of U 4+ . On the other hand, no reflection is ob- 
served at (001) position and also in the scans along the 
principal axes in the first Brillouin zone: (1 + £,0,0), 
(1 + C, 0, 1 - C) and (2, 0, C) for < C < 1- These results 
indicate that the type-I AF structure with moments po- 
larized along the c axis is unchanged by the application 
of a || [100]. The development of the magnetic scattering 
is also observed for a || [110] (Fig.^c)). From the same 
analyses, we confirm that the AF structure is unchanged 
also for cr || [110]. In contrast to the compression along 
the basal plane, the increase of magnetic reflections for 
cr || [001] is very small (Fig.^d)), indicating that the AF 
state strongly depends on the direction of a. 

The magnetic Bragg peaks observed at (100) and (111) 
were fitted by the Lorentzian function convoluted with 
the Gaussian resolution function, to estimate the correla- 
tion lengths £ of the AF moment. The instrumental reso- 
lutions are estimated from the widths (FWHM) of higher- 
order nuclear reflections measured at the corresponding 
Q positions without PG filters. At ambient pressure, £ 
along the [100], [001] and [111] directions are estimated 
to be about 150 A, 260 A and 330 A, respectively. They 
increase rapidly by applying a along the [100] and [110] 
directions. Above 0.3 GPa, the peak widths approach the 
resolution limit (~ 1000 A), and the simple fits give the 
£ values of approximately 2.5 times larger than those for 
cr = 0. On the other hand, £ || [100] for a || [001] remains 
around a small value of ~ 230 A even at 0.46 GPa. These 
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FIG. 1: The uniaxial-stress variations of the magnetic Bragg 
peaks of URu2Si2, obtained from (a) the longitudinal and (b) 
the transverse scans at the (100) position for a || [100], and the 
longitudinal scans at (c) the (111) position for <r|| [110] and 
(d) (100) for o- 1| [001] at 1.5 K. The horizontal bars indicate 
the widths (FWHM) of the resolution limit estimated from 
the higher-order nuclear reflections. Note that the data for 
(7 = are 4 times enlarged in (a), (b) and (c). 



results indicate that the increase of £ is accompanied by 
the enhancement of the AF Bragg-peak intensities. 

Displayed in Fig. |21 is the a dependence of the stag- 
gered moment /x at 1.5 K. The magnitudes of /i G are 
obtained from the integrated intensities of the magnetic 
Bragg peaks at (100) for cr 1 1 [001] and [100], and at (111) 
for a 1 1 [110], which are normalized by the intensities of 
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FIG. 2: Uniaxial-stress dependence of the staggered moment 
fi at 1.5 K. The values of fi are estimated by assuming 
homogeneous AF order. The broken lines are guides to the 
eye. 



the weak nuclear (110) reflection for a\ \ [001] and [110], 
and (101) for a || [100]. We should note that the \i a val- 
ues estimated here are based on the assumption of ho- 
mogeneous AF order. At a = 0, \i is 0.020(4) /is/U, 
which roughly corresponds with the values of previous 
investigations^^ As stress is applied along the [100] 
direction up to 0.25 GPa, \i is strongly enhanced to 
0.22(2) fig/U, and then shows a tendency to saturate 
above 0.25 GPa. The /x value at 0.55 GPa is estimated 
to be 0.25(2) fj, B /U. The fj, a (a) curve for cr 1 1 [100] is 
quite similar to that for the hydrostatic pressureA^ This 
similarity strongly suggests that the enhancement of fi D 
under a is also attributed to the increase of the AF vol- 
ume fraction. However, the estimated rate of increase, 
d/j, /da ~ 1.0 /xs/GPa, is much larger than that for the 
hydrostatic pressure, d/j, /dP ~ 0.25 ^s/GPa. Interest- 
ingly, \i Q also develops with a | [110], tracing the curve for 
a || [100] within the experimental accuracy. For a \ \ [001], 
on the other hand, \i slightly increases to 0.028(3) /xs/U 
at 0.46 GPa, with a small rate d/i /da ~ 0.02 /is/GPa. 

In Fig. 13 we plot the normalized Bragg-peak inten- 
sity 1/7(1.5 K) for cr 1 1 [100] and [110] as a function of 
normalized temperature T/T m , where T m is defined as 
the onset temperature of I{T) as follows. Upon cooling, 
I(T) starts increasing at a temperature T+ and exhibits 
a T-linear dependence below T~ x (< T+). The width, 
8T m = T+ - T~, of this "tail" of I{T) is estimated to be 
2-3 K, and we define T m as the midpoint of T+ and T ? ~. 
Although the experimental errors are somewhat large, 
the cr variations of T m fall in the range of ~ ±1.5 K from 
T m (a = 0) ~ 17.7 K, thereby showing a remarkable con- 
trast with the large a variations of /i . The observed weak 
variations of T m are not inconsistent with the a variation 
of T (dT /da = 1.26 K/GPa), which is obtained from 
the electrical resistivity measurements for a 1 1 [100] For 
a weak stress range a < 0.12 GPa, the I(T) curves for 
both the a directions exhibit unusually slow saturation 



4 



1.2 
1 

0.8 
0.6 
0.4 
0.2 



1 

0.8 
0.6 
0.4 
0.2 





S 

□ k 

o II [100] x ' 
(100) ° 

o (GPa) • 
• 0.12 
□ 0.25 
x 0.37 
a 0.55 



r 

25 

g 20 
h- e 15 
10 

A I 

□ » 



0.2 0.4 0.6 

o (GPa) 



x ft 

f 

o || [110] S ^ 
(111) □ 

a (GPa) 

• 0.10 

p 0.20 

■ 0.38 

a 0.44 

x 0.61 



~i r 

25 

g 20 

H E 15 

x a 10 
□ i 
• * 

□ 



0.2 0.4 0.6 

a (GPa) 



Xa 



0.2 0.4 



0.6 0.8 
T/T 



1.2 1.4 



FIG. 3: Temperature dependence of the normalized Bragg 
peak intensities 1/1(1.5 K) for a\\ [100] (top) and a\\ [110] 
(bottom). Temperature is scaled by T m defined in the text. 
The insets show the a dependence of T m . 



with decreasing temperature. For further compression, 
I(T) shows a sharper onset and more rapid saturation, 
pronounced in a rounding curvature for T/T m < 0.6. 



B. Elastic neutron scattering under 
zero-stress-cooled condition 
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FIG. 4: The uniaxial-stress variations of (a) the staggered 
moment fi and (b) the correlation length £ along the [111] 
direction for <r|| [110] in URu2Si2, measured at 1.4 K after 
cooling the sample at a = 0. The /i data taken under the 
stress-cooled condition (pressurized at 100 K and room tem- 
perature) are also plotted. The inset of (a) shows the a vari- 
ations of [i for a || [100] in U(Ruo.99Rho. 01)2812, measured at 
1.4 K under zero-stress-cooled condition. The lines are guides 
to the eye. 



Figure 0| (a) shows the a variations of ^t for a\ \ [110], 
obtained from the increasing and decreasing a sweeps 
at 1.4 K under the zero-stress-cooled condition. The [i 
value is estimated from the integrated intensity of the 
(111) magnetic Bragg peaks normalized by using the nu- 
clear (110) reflection. The n {p) curve shows a clear hys- 
teresis loop. As a is applied, fi a develops linearly from 
0.016(4) hb/U (a = 0) to 0.20(1) /x B /U (a = 0.4 GPa). 
Upon decompression, on the other hand, \i shows nearly 
cr-independent behavior between 0.4 and 0.3 GPa, and 
then starts decreasing. After the cycle of compression, 
/i Q returns approximately to the initial value at ambient 
pressure. The curve for the cr-decreasing process 

is very similar to that obtained under the stress cooled 
condition. 

In general, the application of a may increase the crys- 
talline mosaic, which weakens the extinction of reflection 
and leads to a significant error in the estimation of the 
intrinsic neutron scattering intensity. Within the pres- 
sure range of the present measurements, the intensity of 



the magnetic (111) peak is always smaller than that of 
the nuclear (110) reference peak. Normally, the stronger 
the reflection, the larger the influence of extinction. We, 
however, observed that the integrated intensity of the 
(110) peak increases by 15%, which is much smaller than 
that of the (111) peak. In addition, the difference of the 
(110) peak intensity between the increasing and decreas- 
ing a sweeps is within the range of 4%, which is also 
much smaller than that of the (111) peak intensity. The 
observed enhancement of the magnetic (111) reflection is 
thus not due to the variation of the extinction effects. 
We also checked the instrumental error of a between 
the stress-increasing and decreasing processes by using 
a strong (220) nuclear reflection. The a variations of the 
integrated intensity due to the extinction effects show no 
significant hysteresis, and we confirm that the error of 
a between the two processes is at most ±0.02 GPa, as 
indicated by error bars in Fig. 0] 

The widths (FWHM) of the (111) magnetic Bragg- 
peaks are slightly larger than the instrumental resolution, 
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and in Fig. 0] (b) we show the correlation length £ of the 
AF moment along the [111] direction as a function of a. 
At ambient pressure, £ is estimated to be about 340 A. 
As a is applied to 0.4 GPa, £ increases to ~ 500 A. Upon 
decompression, it continues to increase, reaches a maxi- 
mum at ~ 0.3 GPa, and then returns to near the initial 
value. Although the experimental errors are large, one 
can see a qualitative correspondence between the \x {<j) 
and the £(<x) curves. The hysteresis loops observed in 
the \i (p) and £(cr) curves strongly suggest that the a- 
induced AF order is metastable and separated from hid- 
den order by a first-order phase transition. 



Strains 



Elastic Constants 



£ = -^{sxx +e yy + e zz ) c al = i(2cn + 2ci2 +4ci 3 +c 33 ) 



yy \ -^12 _ \/2 



(Cll + C12 — C13 — C33) 



e s = y/2e xy 
e\ — \f2e zx 



c a2 = i(cn + C12 - 4ci3 + 2c 3 
c 7 = en - C12 
C S — 2C66 
C 6 = 2C44 



TABLE I: The symmetrized strains and elastic constants in 
the tetragonal symmetry 27 . 



IV. DISCUSSION 

A. Crystal strains under hydrostatic pressure and 
uniaxial stress 

It is important to remark that the uniaxial stress ap- 
plied in the a — a plane brings about similar characteris- 
tics of the AF order, magnitude of \i as well as its T and 
a dependences, to those given by hydrostatic pressured 
This implies that there is an implicit and common pa- 
rameter leading to an equivalent effect in the different 
types of compression. In this subsection, we discuss the 
crystal strains caused by P and cr, and propose that the 
c/a ratio plays an important role in the competition be- 
tween the two ordered phases. 

Within the linear approximation, the uniaxial stresses 
in the tetragonal crystal symmetry are coupled with the 
strains by the relation, 



Let us now consider the influence of the symmetry- 
breaking strains e 7 , e s , e\ and e| on the AF order. It 
is obvious that none of them can be caused by P and 
cr || [001]. On the other hand, a\\ [100] and a || [110] give 
rise to e 1 and e s , respectively. Therefore, if the AF order 
is induced by the symmetry-breaking strains, it should 
occur only by a || [100] and a || [110], and it is not neces- 
sary for their effects to be the same. This is inconsistent 
with our experimental results: /i is induced by both 
P and a (in-plane), and a || [100] and a || [110] have the 
same effects within the experimental accuracy. We thus 
conclude that the symmetry-breaking strains are irrele- 
vant to the evolution of the AF phase, at least, in the 
weak pressure range. 

We next consider the variations of the symmetry- 
invariant strains, e al and e a2 , which can be expressed 
by the relative variations of the unit cell volume, V, and 
the c/a ratio, 77, as follows: 
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where the ct^'s, Cy's and e^-'s indicate the uniaxial 
stresses, elastic constants and strains. The elastic en- 
ergy symmetrized in the tetragonal point group can be 
expressed by the formi^i 

E el = l c al (s a1 ) 2 + c al2 e al e a2 + ^c a2 (e a2 ) 2 (2) 

+ ^ 7 (e 7 ) 2 + icV) 2 + \c e {{e\f + (el) 2 }. 

The definition of the c*'s and e"s are given in Tabled 
These notations for the strains are useful in discussing the 
symmetry of lattice distortion. For example, the strains 
of e 7 , e 5 , and e e types break the tetragonal symmetry, 
while the strains of e a type change the volume and the 
c/a ratio, conserving the tetragonal symmetry. We also 
show in Table [H] the symmetrized strains divided by the 
stresses, e'/X for X — P, a\\ [100], a \\ [110] and a \\ [001], 
calculated from Eq. QJ. 



t t- &XX "!" Syy + S ZZ A/36 , (3) 

Vq 

v - m 



— ^22 



Vn 



(4) 



where Vq and rjo denote the values at ambient pressure. 



Using the known ctj values of URu2Si2 (Table HID, we 
calculated the rates of change in the volume, dv/dX, 
and the c/a ratio, dfj/dX, in Table IIVI The calcula- 
tions show that 77 is increased by a\\ [100] and [110] at 
the same rate, dfj/da ~ 3.0 x 10 -3 GPa -1 . Interest- 
ingly, fj is also expected to increase under hydrostatic 
pressure, because of the Poisson's effect. From the cal- 
culations we obtained the relation between the increas- 
ing rates: dfj/da ~ 3 x dfj/dP. These features seem to 
be consistent with the experimental results that \x are 
equally enhanced by a || [100] and [110], and the relation 
djio/da ~ 4 x dfi /dP holds. The observed /i D (-P) and 
fio(c) curves are well scaled by 77 (Fig.EJ, indicating that 
the c/a ratio is relevant to the competition between the 
two types of order. On the other hand, the volume con- 
traction v is irrelevant, because P should exert a stronger 
influence than cr, which is inconsistent with the obser- 
vation. In this context, however, fi a is expected to be 
suppressed by applying er|| [001], whereas actually it is 
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almost independent of the stress (Fig. [2J. This can be 
understood, if the AF phase observed at ambient pres- 
sure is caused by irremovable local distortions which are 
"pinned" near impurities and defects. 
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FIG. 6: The uniaxial-stress (a || [100]) variations of the stag- 
gered moment fi for U(Ruo.99Rho.oi)2Si2, measured at 1.4 K 
after cooling the sample at a = 0. 



FIG. 5: The spatially averaged AF moment fi obtained from 
the elastic neutron scattering under hydrostatic pressure P 
and the uniaxial stresses a || [100] and [110], plotted as a func- 
tion of fj = (r] - T] )/r}o = e zz - (l/2)(e xx + e yy ). 

The magneto-elastic energy E me for the type-I AF or- 
der in the tetragonal crystal is given by 

E me = -D v vM 2 - D v fiM 2 , (5) 

where M denotes the staggered magnetization and D v ^ 
magneto-elastic coupling constants^ The above consid- 
eration implies that \D V \ is larger than \D V \ in URu2Si2- 
This is supported by recent thermal-expansion measure- 
ments performed under P, which revealed that the c/a 
ratio significantly increases as the AF phase develops 
with decreasing temperaturei^S* 3 ^ 

The significance of the c/a ratio is also recognized from 
the behavior of the alloy system U(Rui_ :c Rh. E )2Si2. In 
this system, the c/a ratio is known to increase as x 
increases^? For x ~ 0.02, fj reaches ~ 1 x 10~ 3 : the 
value at which the AF phase is fully induced in the pure 
compound (see Fig-JSJ. Correspondingly, the AF phase is 
found to develop at x ~ 0.015. 33 To test the relevance of 
the "chemical stress" to the phenomena, we applied uni- 
axial stress (|| [100]) to the alloy U(Ruo.9gRh .oi)2Si2. 
We observed that /J (T = 1.4 K) steeply increases with 
a, from 0.026(3) to 0.20(2) /jb/U, and the saturation of 
fi is more abrupt than that for the pure system (Fig.^JJl. 
These facts indicate that the axial strain, which is gener- 
ated by Rh doping, also governs the two phase competi- 
tion in this alloy system: the Rh 1% system is chemically 
compressed near to the AF instability point, already at 
ambient pressure. The hysteretic behavior is also de- 
tected in the /x (cr) curve, supporting the argument that 
the transition is of first-order. 

In our previous measurements using hydrostatic pres- 
sure, we observed a sudden increase in [i Q from ~ 0.22 to 
~ 0.40 /xs/U at P c ~ 1.5 GPa. If this anomaly is also 



caused by the increase in fj, then similar behavior should 
be observed at a (_L [001]) ~ 0.6 GPa, where fj is expected 
to reach the value (~ 1.8 x 10~ 3 ) estimated at P c . The 
maximum applied a in the present study is 0.61 GPa (|| 
[110]), and in this a range we observed no indication of 
the P-transition (see Fig. Eland Fig. El upper right data 
points). Further investigation with higher stress will be 
needed to resolve the origin of this anomaly. 



B. The application of the Landau theory 

The stress-induced first-order phase transition ob- 
served in URu2Si2 is qualitatively understood in terms 
of the Landau's free energy theory with a time-reversal- 
invariant order parameter as follows. We assume the free 
energy F(ip, M) of the form, 20 i 31 i 34 i 35 



F = -r i> i) 2 + u^ i + -r M M 2 + u M M 4 (6) 

+ 2m v ,mV' 2 M 2 , 
= Gty,(T - T^), (7) 
ru = a M (T -T M ), (8) 

where ip and M denote the hidden order parameter and 
the staggered magnetization, and the signs of at and Ui 
are positive. It is straightforwardly seen that a first or- 
der phase transition between ip and M may occur at 
the boundary r$ — \Ju^/umTm (< 0) on the condi- 
tion u^um < iQj M . Suppose that only the symmetry 
invariant strains are relevant. Then the total free energy, 
F t = F + F e i + F me , including the elastic energy F e i and 
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X 
Strain 

e al /X 

e a2 /X 
e"</X 

£ S /X 

ej/X,e\/X 



{-P,-P,-P, 0,0,0) 

\/3 -2cJ 3 + (cii+c 12 )c33 

2 eil+c 12 — c 13 — C33 

3 - 2c 13 + ( c ll + c 12)c33 







g || [100] 
(-a, 0,0, 0,0,0) 

1 — C13+C33 

v' 3 -2cJ 3 +(c 1 i+c 12 )c33 

1 2C13+C33 

VO -2cj 3 +(c 1 i+c 12 )c33 

1 1 

V2 cu-c 12 





gj [110] 

(-a/2, -a/2, 0,0,0, -a/2) 

1 -C13+C33 

\/3 - 2c i 3 + ( c ll + c 12)c33 

_L 2C13+C33 

\7S - 2c 13 + ( c ll + c 12)c33 


1 1_ 





g II [001] 
(0,0, -a, 0,0,0) 

__1 C11+C12 — 2ci3 

v / 3_-2cJ 3 + (c 1 i+c 12 )c33 

2 C11+C12+C13 

3 - 2c 13 + ( c ll+ c 12)c33 







TABLE II: The symmetrized strains divided by stresses induced by hydrostatic pressure and uniaxial stresses in the tetragonal 
symmetry. 



Cll 


C33 


C44 C66 


Cl2 C13 


25.5 


31.3 


13.3 18.8 


4.8 (8.6) 






(xlO 11 


erg/cm 3 ) 



TABLE III: The elastic constants at low-temperatures 
for URu2Si2 obtained by the ultrasonic-sound-velocity 
measurements^- The value for C13 was estimated from a com- 
parison between URu2Si2 and the isostructural compounds 
CeCu2Si2 and CeRu2Si2- 



X 


P a 


1 [100] 


a|| [110] a|| [001] 


dv/dX 


-7.3 


-2.8 


-2.8 -1.6 


df]/dX 


1.2 


3.0 


3.0 -4.9 








(xlO -3 GPa" 1 ) 



TABLE IV: The increasing rate of the symmetry invariant 
strains induced by hydrostatic pressure and uniaxial stresses, 
calculated from the elastic constants of URu2Si2- 



the magneto-elastic energy F me becomes, 

F t = ^r^ 2 + u tf ,ip i + ^r / M M 2 +u M M 4 

2 t\,t2 



6 3 3 



(9) 



r' M = a M (T (T M + ^v + ^ry)). (10) 

Here, we neglected the coupling between ip an d the 
strains for simplicity, but it should be remembered that 
tp seems to be also weakly coupled to rj. This is ex- 
pected because r\ increases below T , 36 and because T a 
increases with P 37 i 38 i 39 i 40 i 41 and a(\\ [100])^ T$ must be 
larger than Tm at ambient pressure, since hidden order 
forms the majority phase. If D n > and ID^I 3> \D V \ 
in this situation, then the AF transition temperature T' M 
(= T«+— w+— n) increases with increasing fi, so that 
the first order phase transition occurs at the critical point 
f] c as shown in Fig. EI a). By comparing the expected 





Para 


(a) 


2nd 








1st 


\|» = 


M = 




M*0 





N(f|) P= a = 



A 

11c 




FIG. 7: (a) The schematic drawing of the fj — T phase diagram 
expected from the Landau's free energy consideration involv- 
ing the elastic and magneto-elastic interactions. The phase 
boundary between hidden order and the AF order (7) = is 
characterized by the first order phase transition. The f) c value 
of URU2S12 is estimated to be ~ 10 -3 . (b) The schematic 
drawing of the distribution of 77. The width of distribution is 
expected to be ~ 10 -4 . The distribution is shifted right by 
thermal expansion caused by hidden order (the broken line) 
and by compression (the dotted line) , generating the inhomo- 
geneous AF phase. 



phase diagram with the present experimental results, r) c 
is roughly estimated to be ~ 1CP 3 in URu2Si2. Since 77 
is an extensive variable, in principle the phase diagram 
should have an area near f] c where r) shows a discontin- 
uous change. Such an area is, however, expected to be 
very narrowj 3 ^ 3 ^ and not described in Fig. 7. 

The above consideration is intended for a homogeneous 
system, and does not account for the inhomogeneous de- 
velopment of the AF phase. The crucial feature would 
be the smallness of f) c . Here we suggest the presence of 
random distribution of r] in URu2Si2, due to some im- 
perfection of the crystal, as schematically shown in Fig. 



8 



Efb). The width of the distribution is expected to be 
of the order of 10~ 4 , which will be hard to detect and 
analyze using the usual microscopic probes. At ambient 
pressure, the mean value of fj (taken as in Fig. Efb)) 
should be smaller than fj c , so that most part of the sam- 
ple shows hidden order below T^. We should remember 
here that the linear thermal-expansion coefficients show 
an increase of fj of the order of 10 -4 below T^£& fj is 
thus expected to exceed fj c in small fragmentary regions 
of the sample, where the AF order takes place, being 
detected as tiny moment on volume-average. By apply- 
ing P or <7_L [001], the mean value of fj exceeds fj c , and 
the AF volume fraction inhomogeneously develops to the 
whole part of the sample, as is observed in the 29 Si-NMR 
measurements under P.^ilfl The temperature and stress 
dependence of the AF volume fraction should strongly 
depend on the condition of sample preparation, because 
such has a strong influence on the compressibility, the 
thermal expansion, and the distribution function of fj. 
This is consistent with the observed annealing effects, 
where the magnitude, the onset temperature and the T 
variation of the AF Bragg-peak intensity all show signif- 
icant sample-quality dependence^ In particular, in this 
context the onset temperature of I(T), which we define 
as T+ in this paper, could become higher than T a , if the 
distribution of fj extends over fj c above T a . This is actu- 
ally observed in the present system^MiLiSi^ where the 
width of onset |T+ — T~| strongly depends on the specific 
experiment and sample. We emphasize that the AF re- 
sponse of such variety of starting conditions at ambient 
pressure is dominated by undetectably small change in 
the c/a ratio. 

Through the above considerations, we have stressed 
that the weak magnetism at ambient (and very low) pres- 
sure could reasonably be understood as the mixing of 
the high-pressure AF phase. This allows ones to adopt 
a time-reversal-invariant hidden-order parameter such as 
quadrupole moment. However, the presence of the AF 
fraction at very low pressure has not yet been confirmed 
by experiments. The present experiments do not ex- 
clude the possibility that the low-pressure magnetism is 
induced by an order parameter that breaks time reversal 
invariance but is nearly non-magnetic, such as an oc- 
tupole moment^ 



V. CONCLUSION 

We have presented elastic neutron scattering experi- 
ments under uniaxial stress on single crystal URu2Si2, 
and discussed the nature of the unusual competition be- 
tween hidden order and inhomogeneous AF order. A 
significant increase of the AF Bragg-peak intensity was 
observed when a is applied along the [100] and [110] axes, 
while it is nearly constant for <r || [001]. The a varia- 
tion of the AF scattering intensity for a\ \ [100] roughly 
corresponds with that for a\ \ [110], indicating that the 
AF evolution is isotropic with respect to compression 
in the tetragonal basal plane. The isothermal curve of 
the AF Bragg-peak intensity, which was obtained for 
U(Ruo.9gRho.oi)2Si2 as well as URu2Si2 under the zero- 
stress-cooled condition, shows a clear hysteresis loop, in- 
dicating that the phase transition from hidden order to 
the AF order is of a first-order. It was also found that 
the application of uniaxial stress enlarges the AF phase 
more effectively than that of hydrostatic pressure. We 
considered the crystal distortions induced under a and 
P, and pointed out that the observed features can rea- 
sonably be explained by the increase of the c/a ratio with 
the compression. This interpretation is consistent with 
the results of the recent thermal-expansion measurements 
performed under hydrostatic pressure.^ The inhomoge- 
neous development of the AF phase can also be ascribed 
to the presence of random axial strains with a very small 
distribution width of ~ 10~ 4 . 
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